A reality proof in "PT-symmetric quantum mechanics 
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We review the proof of a conjecture concerning the reality of the spectra of certain "PT- 
symmetric quantum mechanical systems, obtained via a connection between the theories 
of ordinary differential equations and integrable models. Spectral equivalences inspired by 
the correspondence are also discussed. 
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1 Introduction 

The unexpected reality of the spectra of a class of non-Hermitian quantum 
mechanical systems has provoked a fair amount of work since the first example 
was noted by Bessis and Zinn- Justin ,1 . Generalising the result at M = 3/2, a 
numerical study of the Schrodinger equations 

led Bender and Boettcher to conjecture that under suitable boundary conditions 
the spectrum of 0) was entirely real and positive provided M> 1 [|. 

Bender and Boettcher also noticed the relevance of 'PT-symmetry to these 
problems. For all M, the potentials V(x) — —(ix) 2M are invariant under the com- 
bined action of parity and time reversal, and this implies the eigenvalues are either 
real or occur in complex-conjugate pairs. If all energy eigenstates are also eigen- 
states of the VT operator, then the spectrum is guaranteed to be entirely real, 
but this is not true in general: the PT-symmetry can be spontaneously broken, in 
which case the corresponding energies occur in complex-conjugate pairs. In fact, a 
proof of the reality conjecture has only recently been obtained [21, via a correspon- 
dence between ordinary differential equations and integrable models (the so-called 
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Fig. 1. The initial 'phase diagram' at fixed M. 



'ODE/IM correspondence' (aEHZIIal)- Here we review the proof and touch on a few 
of the byproducts of the correspondence relevant to the ODE community. 

Generalising the Bessis - Zinn- Justin and Bender - Boettcher Hamiltonians, we 
consider the class of 'PT-symmetric problems 



dx 2 



(%xf M -a(ix) M - L + 



1(1 + 1) 



il>{x) = X^j(x) , 1>(x) G L 2 (C) . (2) 



For M < 2, the contour C may be taken to be the real axis with a small dip below 
the origin if 1(1 + 1) ^ 0. For larger M, the contour should be distorted into the 
complex plane to ensure the correct analytical continuation of the original problem, 
as explained in [2] . The result of 0] is that the spectrum of is real for M > 1 
and a < M+ 1 + \2l+ 1| and positive for M > 1 and a<M+l-|2Z+l|. Referring 
to figure n reality was proven for (a, I) 6 B U CU D, and positivity for (a, I) € D. 



2 The functional relation 



Rather than considering the PT-symmetric problem alone, we introduce a 
related problem obtained from J2J) by sending x — * x/i and E — > — E: 

(-^ + x™ ± ax^ + ^) *(*) = • (3) 

It will be convenient to treat positive and negative values of a together. The bound- 
ary conditions for are that should vanish as x — > oo along the real axis, 
and behave as x l+1 as x — > 0. For Z real and larger than —1/2, this problem is Her- 
mitian. In the language of ordinary differential equations in the complex domain, it 
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is sometimes called a 'radial' problem, while is called a 'lateral' problem. (Note, 
the radial problem can also be considered for $lel < —1/2, but is best then defined 
by analytic continuation in I.) 

Following the ideas particularly advocated by Sibuya and Voros in the con- 
text of ordinary differential equations [51ll(J| , we approach the spectral problems by 
studying the behaviour of the associated spectral determinants. Adopting the con- 
vention used in 15 , let {E^'} be the set of eigenvalues of with inhomogeneous 

term =Fa(ix) M ~ 1 , and let {e^} be the eigenvalues of |JS} with inhomogeneous term 
±ax . Then define two pairs of spectral determinants, as follows: 

r (±) ( £ )=T W(0)n (i + JLj) . w 
j=o \ E j I 



and 



Q^(E)=Q^(0)f[[l--^) . (5) 



Both products are convergent for M > 1, and define entire functions of E, and 
their zeroes (or, for T(E), their negatives) coincide with the eigenvalues of the 
corresponding spectral problem. For M < 1 convergence factors must be added, 
and it is more efficient to define the spectral determinants indirectly, via certain 
special 'Sibuya' solutions to (J3J, which are anyway needed to prove the key identity 
(jHJ below; see [HIE] for more details. 

By considering the asymptotic behaviour of the Sibuya solutions in the complex 
plane it is not hard, following the arguments given in 0, to establish a Stokes 
relation, from which one can obtain the following functional equations 

T^\E)Q^(E) =Lo-^ l+l±a ^ 2 Q^{u J - 2 E)+J 2l+l±a ^ 2 Q^\u J 2 E) , (6) 

where to = e l7r /( M + 1 ). This shows that the spectral problems J2J and J21 are related 
by much more than a simple change of variables and boundary conditions, and 
also that the spectral problem at positive a is necessarily tied up with the problem 
at negative a. Furthermore, this functional relation taken at a = is well-known 
in the integrablc model (1M) world, where it goes by the name of Baxter's TQ 
relation. This observation, initially developed in |5l(il7) . has led to an exact mapping 
between ODE quantities such as spectral determinants and functions constructed 
in the context of integrable models. It provides a powerful tool to analyse ODE 
problems using IM techniques, and vice-versa. 

If we set E = in J^J we can use the entirety of T and the product form for 
Q © to obtain 

oo / o (T)_,„2 (±) \ 

(7) 

In IM language sets of coupled equations of this type are known as Bethe ansatz 
equations. Their solution allows Q to be reconstructed via © and thus T using 




Czech. J. Phys. 53 (2003) 



A 



3 



Dorey et ad. 



the TQ relation © . Finally the eigenvalues of the 'PT-symmetric Hamiltonian are 
found by searching for the zeros of T{E). Alternatively, a technique developed in 
integrable models allows one to convert the Bethe ansatz equations J7J into a single 
nonlinear integral equation, from which one can easily numerically obtain both the 
zeros of Q(E) and T(E), as first emphasised in [5J[7|. 



3 The proof 

Returning to the reality proof, a second set of Bethe ansatz like equations can be 
obtained from @ if we instead set E = — Ey~' and invoke the entirety of Q^~\E). 
Rearranging and using the product form floj leads to 



e 



E 



This equation couples the so-far mysterious eigenvalue Ej of the PT-symmetric 
problem J2J) to the much better-controlled eigenvalues of the Hermitian problem 
(|2J). Indeed, a Langer transformation JT] shows that the eigenvalues e^ - ) of the 
radial problem are also positive if a < M + 21 + 2 ^ . It is from this equation that 
we are able to prove the reality of the PT-symmetric problem. 

The first step will be to set E^ — \Ej\ exp(iSk)- Now take the modulus 2 of 
(JSJ to obtain 

fr / (^ ) ) 2 + l^ +) | 2 + 2^ ) |^ +) |cos( I |f T +^-) \ 
111 (ej- J )» + |J5< + >P + 2 eW|^ + >l-' ^ rJ ' 



■ | ^wo V M+1 



For a < M+ 21 + 2 all the e\ are positive, and each single term in the product on 
the LHS of @ is either greater than, smaller than, or equal to one depending only 
on the relative values of the cosine terms in the numerator and denominator. These 



are independent of the index i. Therefore the only possibility to match the RHS is 
for each 
requires 



for each term in the product to be individually equal to one, which for E [ + ] ^ 



™<WTi + S i) = C0S (STI - *i) . OT sin ( aTTi) sin (^') = • (!0) 
Since M > 1, this latter condition implies 

5j = mr , n£Z (11) 

and this establishes the reality of the eigenvalues of @ for M > 1 and a < M+21+2 
or, relaxing the condition on I, a < M + 1 + |2Z+1|. 

At M < 1 the simple form JSJ for D^~'(E) acquires extra convergence factors, 
resulting in a break down of the proof, as expected given the numerical findings 
of 0C1 More details can be found in 0j; for a further generalisation of the method, 
see H2). 
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Fig. 2. The domain of unreality for M = 3. Horizontal and vertical axes are I and a 

respectively. 

Finally, we remark that while the above constraints on the parameters M, a 
and I are sufficient conditions they are not necessary, as demonstrated in figure [2] 
for the case M = 3. The full domain of unreality, obtained numerically in |13j. 
is the interior of the curved line, a proper subset of A which only touches its 
boundary at isolated points. In the small, approximately-triangular region inside 
A but outside the curved line which abuts the crossing of the lines separating the 
regions A, B, C, D, the spectrum is not only real but also entirely positive, despite 
the fact that it lies outside the domain D. The separating lines themselves have 
further significance, in that at these points the model has a hidden supersymmetry, 
which plays a significant role in the breakdown of reality |13j . 

4 Further consequences 

One further application of the ODE/IM correspondence was also found in 4 , 
and concerns the lateral and radial problems at M = 3. As remarked in it is 
convenient to reparametrise the angular-momentum term by setting A = y/3(2l+l) , 
so that the equation becomes 

(-5? +x6+ ax2 + w) * {x) = mx) • (12) 

Define a = (a, A) T and let D(E, a.) denote the spectral determinant for this prob- 
lem. Via the Bethe ansatz approach, it turns out that this problem has a relationship 
with a third-order ordinary differential equation: 

(D( ff2 -2)D( 5l -l)D(.g ) +3»)4>=^.E<I> (13) 
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where D(g) = (£ - f ), and 

,g = 1 + (a + V3A)/4 , g x = 1 + a/2 , g 2 = 1 + (a - \/3A)/4 . (14) 

This third-order equation is associated with S[/(3) Bethe ansatz equations, as dis- 
cussed in |14I15| . Furthermore, the third-order equation is symmetrical in {170, Si, 32}, 
a feature which is completely hidden in the original second-order equation. By play- 
ing with this symmetry, one can establish some novel spectral equivalences between 
different (second-order) radial problems, and also between these and certain lateral 
problems. We add to the results discussed in 0] the remark that, when expressed 
in terms of the parameters (a, A), the mappings turn out to act as certain 2x2 
matrices in the Weyl group of SU(3). The matrices L and T defined as 

generate the Weyl group of SU(3), the matrix L describing a rotation and T a 
reflection. In 0] we found a spectral equivalence between a pair of radial problems 

j(a)D(E, a) = j{Ta)D(E, Ta) , (16) 

where 7(a) = 2v / i7r/r(A/2v / 3). An equivalence was also obtained at a particular 
set of points for which the radial problem D(E, a) has a quasi-exactly solvable 
sector [TD1IT7] . For positive integer J and aj = —(4 J + A/ \/3, A), the first J levels 
of D(E, a.j) can be computed exactly as the zeros of the Bender-Dunne polynomials 
Pj{E) [HI. The result is that modulo the QES levels, the problems D(E, TLcxj) 
and D(E,atj) are isospectral. 

A further spectral equivalence particularly relevant to the main theme of this 
paper occurs between a radial problem and the related lateral problem J5J) taken 
at M = 3. The equivalence is 

T(-E, a) = j(LTa)D(E, LTa) . (17) 

If we denote LTa = (a, A) T then the problem on the RHS is Hermitian for A > 0, 
and this relation provides a simple explanation for the reality of the spectrum of 
the "PT-symmetric problems in these particular cases. At the points a = olj we can 
combine the dualities to prove a simple result concerning the spectrum of T(E, ocj), 
namely that the only energies to become complex as the parameters (a, I) move into 
region A through QES values lie in the solvable part of the spectrum ^1 • This adds 
to a similar, as yet unproven, conjecture concerning quartic QES potentials |22j . 

More detailed reviews of the ODE/IM correspondence, with more extensive 
sets of references, can be found in |19U2l)ll2*T] . The original correspondence relates 
only the groundstate of a quantum field theory to an ODE, but recently a set of 
differential equations has been found for the excited states of the field theory |23j . 
Further aspects of the correspondence are still being developed, and much remains 
to be understood. In particular, it would be valuable to have a more physical insight 
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into why the relationship between integrable quantum field theories and ordinary 
differential equations should be so close. As more examples are uncovered, we can 
start to hope that progress on this so-far mysterious issue may not be too far away. 
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